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Abstract 

We consider the Lame system of linear elasticity when the inclusion has the ex- 
treme elastic constants. We show that the solutions to the Lame system converge in 
appropriate H -norms when the shear modulus tends to infinity (the other modulus, 
the compressional modulus is fixed), and when the bulk modulus and the shear mod- 
ulus tend to zero. Using this result, we show that the asymptotic expansion of the 
displacement vector in the presence of small inclusion is uniform with respect to Lame 
parameters. 
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1 Introduction 

Recently there is growing interest in partial differential equations with high contrast coeffi- 
cients in various contexts. Among them are the photonic and phononic band gap problems 
where the electromagnetic parameter or the bulk modulus tend to infinity, biomedical imag- 
ing where anomalous tissues have large material parameters, and the stress concentration in 
between two inclusions with extreme material properties to name a few. See [10l 13] and 
references therein. The purpose of this paper is to prove two basic theorems in relation to 
the PDEs with high contrast coefficients. The first one is to show that when the material 
parameters tend to the extreme, the corresponding solutions converge strongly in appropri- 
ate norms. The other one is to show that the asymptotic expansion of the solution in the 
presence of small inclusions holds uniformly with respect to material parameters. We prove 
these facts in the context of the system of linear elasticity. Corresponding results for the 
conductivity equation (a scalar equation) have been obtained in [32J H31 US] ■ 

We consider a linear isotropic elastic body containing an inclusion with different elastic 
parameters. When the bulk and shear moduli of the inclusion are finite, the solution satisfies 
the transmission condition along the interface (the boundary of the inclusion) . If the shear 
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modulus of the inclusion is infinity, then the interface transmission condition is replaced by 
a null condition of the displacement (see Section 2). If the bulk and shear moduli are zero, 
then it is replaced by the traction zero condition on the boundary of the inclusion. The 
first objective of this paper is to prove the convergence in an appropriate H 1 space of the 
solution to the Lame system as the bulk and shear moduli tend to the extreme (zero or 
infinity) (see Theorem 14. 2 1) . 

The second objective of this paper is to prove a closely related problem of uniformity of 
the asymptotic expansion. In imaging small inclusions from boundary measurements, it is 
of fundamental importance to catch the boundary signature of the presence of anomalies. In 
this respect, an asymptotic expansion of the boundary perturbations of the solutions due to 
the presence of the inclusion, as the diameter of the inclusion tends to zero, has been derived. 
The asymptotic expansion is derived in [22j [18j |6] for the conductivity (scalar) equation, 
and in [T2l [28] for the Lame system of linear isotropic elasticity. The asymptotic expansions 
have been effectively used for imaging diametrically small inclusions. See for example [181 
|2"51 [71 [HI [23 E] • We also mention the topological derivative based shape optimization where 
the asymptotic expansion is an essential ingredient (see for example [TTJ [531 12] ) • hi [2] 
topological derivative based detection algorithms for the localization of an elastic inclusion 
of vanishing characteristic size have been developed and their resolution and stability with 
respect to measurement and medium noises analyzed. 

In these applications, it is important to know that the asymptotic expansion holds uni- 
formly with respect to the pair of Lame parameters. We prove this in the second half of this 
paper under the assumption that the compressional modulus is bounded, which is necessary 
(see Theorem 16.11 for precise statements). It is worthwhile to mention that this result may 
have a relation with the cloaking as discussed in 29. 114]. 

The methods of this paper are different from those of [29], where uniform validity of 
the asymptotic expansion for the conductivity (scalar) equations is proved, in that they are 
based on the layer potential techniques. The solutions to the Lame system can be expressed 
as a single layer potential on the boundary of inclusion. We show that if -1 / 2 -norms of the 
potentials are bounded uniformly with respect to Lame parameters, and the main results 
follow from this fact. 

This paper is organized as follows. In section 2, we set up the problems for finite and 
extreme moduli, and review the representation of solutions using layer potential techniques. 
In section 3, we prove that the energy functional is uniformly bounded. As a consequence, we 
obtain that the potentials on the boundary of the inclusion are uniformly bounded. In section 
4 we show that these potentials converge as the bulk and shear moduli tend to extreme 
values and prove Theorem 14. 2 1 In section 5, we briefly discuss that similar boundedness and 
convergence result hold to be true for the boundary value problem. Section 6 is to prove 
Theorem 16.11 which asserts that the small volume expansion holds independently of Lame 
parameters. The results and methods hold to be true even if there are multiple inclusions. 
We make a brief remark on this in the last section. 

2 Problem setting and representation of solutions 

Let D be an elastic inclusion which is a bounded domain in M. d (d = 2, 3) with the Lipschitz 
boundary. Let (A, /i) be the pair of Lame (shear and compressional) parameters of D while 
(Aq, no) is that of the background R d \ D. Then the elasticity tensors for the inclusions and 



2 



the background can be written respectively as C 1 = iC\j ks ) and C° = (C°- M ) where 
and the elasticity tensor for M. d in the presence of the inclusion D is given by 



1 C X + (1-1 



D 



(2.1) 



where Id is the indicator function of D. We assume that the strong convexity condition 
holds, i.e., 

jtx > 0, dX + 2[i > 0, fi > 0, and dX + 2^ > . (2.2) 
We also assume that 

(A-A )(/i-/io)>0 ! (2.3) 

which is required to have the representation of the displacement vectors by the single layer 
potential in the following. We also denote the bulk modulus by re which is given by re = 
X + 2(j,/d. 

We consider the problem of the Lame system of the linear elasticity: For a given function 
h satisfying V • C°V s h = in R d , 



f V • CV s u = in R d , 

\ u(sc) - h(x) = 0{\x\ 1 ~ d ) as |x| ->•«>, 

where V s u is the symmetric gradient (or the strain tensor), i.e., 

V s u := ~(Vu + (Vu) T ) (T for transpose). 



(2.4) 



Let 



Cx.^u := V • C x V 5 u = ^Au + (A + ,u)VV • u 



and define the corresponding conormal derivative du/dv on dD by 

^ := C'fV'u) n = A(V • u)n + /i(Vu + (Vu) T )n on dD, 
ov 



(2.5) 



where n is the outward unit normal to dD. Let D c := R d \ D. Let C\ 0tfl0 and be those 
corresponding to (Aq,/Xo)- Then (|2.4[) is equivalent to the following problem: 



'£a ,mo u = ° in -D c , 
Cx.^u = in D, 
u|- = u|+ on dD, 



du 



dv 



du 

dv Q 



(2.6) 



on dD, 



k u(sc) - h(x) = Odccl 1 d ) as |sc| -> oo 



where the subscripts + and — indicate the limits from outside and inside D, respectively. 

We also consider the two limiting cases of when both re and \i tend to 0, and when 

fi — > oo while A is fixed. In relation to the latter case it is worth mentioning that if A — > oo 



3 



and /i is fixed, (|2.6[) approaches to a different problem. Roughly speaking, if A — > oo, then 
V • u is approaching to while AV • u stays bounded. So (|2.6I) approaches to the modified 
Stokes' problem /iAu + Vp = with p = AV • u. (See [J].) Hence we assume that A is 
bounded throughout this paper. 

If k = fi = (or A = /i = 0), one can easily see what the limiting problem should be. 



Since f£H- = 0, we have from the fourth line of 
equation in this case is 

£\ o ,h > vl = ' m D 
du 



that 



du 1 + 



dv 
u{x) 



= 
- h(x) 



on dD, 
= 0(\x\ 1 -*) 



0. So the elasticity 



(2.7) 



as \x\ — )• oo. 



To describe the equation when fi — oo while A remains bounded, we need to introduce 
the following functional space: Let "J be the d(d + l)/2 dimensional vector space defined by 

f:={i/; = (i/; (1) ,--- > W ) T : d^ + drf® = 0, 1 < i, j < d }. (2.8) 

We emphasize that ^ is the space of solutions £a, m u = in D and du/dv = on dD for 
any (A, fi). Let ipj, j = 1, . . . , d(d + l)/2, be a basis of ^. If fi — > oo, then from the second 
and fourth equations in (|2.6[) we have 

Au + VV-u = in D, (Vu + (Vu) T ) n = on dD, 



which is another elasticity equation (with \i = 1 and A 
boundary. Thus there are constants <Xj such that 



0) with zero traction on the 



w(x) 



d(d+l)/2 

E 



ajipj(x), x e D. 



So, the elasticity problem when /i = oo is 



^Ao,mo u = ° in -D c . 

d(d+l)/2 

u = a j^Pj on dD, 



(2.9) 



k u(x)-h(x) =0(|a;| 1 - d ) as \x\ -+ oo. 

Note that the solution u to 



We need extra conditions to determine the coefficients a 
satisfies 



du 

,,d 



■ tpl da — 0, 1 = 1, 



d(d+l) 



(2.10) 



So, by taking a (formal) limit, one can expect that the solution u to (|2.9p should satisfy the 
same condition, and the constants ctj in (|2.9p are determined by this orthogonality condition. 

We now review the representation of the solution to (|2.6p using the single layer potential 
for the Lame system following [27, 20 [5TJ [!|] . The Kelvin matrix of the fundamental solution 
r = (Tij)fj =1 to the Lame system is given by 

P XlX ' if d = 2 , 

x ^ 

if d = 3 , 



r«(aj) := 



2tt 



-<5 r , In Ice I 



47T Ixl 



2tt |a;| 2 

p X-iXj 



2tt |a;l 3 
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where 

1 f 1 1 \ J a 1 f 1 

a = — I — h and p = — I 



2 \/i 2/i + A/ 2 \/i 2fj, + \ / 

When D is a simply connected domain, the single layer potentials of the density function 
tp = Upi ■ ■ ■ (fid] on 3D associated with the Lame parameters (A, //) are defined by 

S D [ip](x):= f T{x-y)<p{y)da{y) , xeR d . (2.11) 

J 3D 

We denote by r°, <5>£>[</?] the fundamental solution and the single layer potential associate 
with the Lame parameter (Ao, ^o) respectively. The conormal derivative of 5d[v] enjoys the 
jump relation on 3D: 

■^Sd[<p]\ ={±^I + K* D )[if] a.e. on9D, (2.12) 

where JC* D is defined by 

f d 

K* D [<p](x) = -p.v. —T(x-y)tp(y)da(y) a.e. x € 3D, (2.13) 

J 3D ° v k 

where p. v. stands for the Cauchy principal value. We denote by K,^ the operator corre- 
sponding to (Ao,/io)- 

We introduce a weighted norm, ||u||#i (q), in two dimensions: let fi be either M. d or D c , 
and let 

/ |ffi!L to+ / \Vu(x)\ 2 dx. 
Jn y/1 + \x\ 2 Jn 

This weighted norm is introduced because the solutions u satisfies only u(x) = 0(|x| _1 ) 
in two dimensions as \x\ — > oo. For convenience in presenting results of this paper, we put 
VK(f2) := H^(fi) in two dimensions, and W(f2) := i? 1 (fi), the usual Sobolev space, in three 
dimensions. Let \fr be the space introduced in (|2.8p . and define 



u llffi,(fi) : - 



H~ 1/2 (dD) :={tpe H- 1 l 2 {dD) d : {ip,^>) = for all V € * }. (2.14) 

Here (<p, tp) denotes the U^^-H 1 ! 2 product. Then ±\l + K* D is invertible on H^ 1/2 (dD). 
We also have 

l|S£>Mllw(R d ) < C\\ip\\ H -i,2 [dD) (2.15) 

for all <p <= H- l / 2 {dD) d . 

It is proved in pi] that the solution u to (|2.6|) is represented as 



h(,)+ 5 o M(x)) x€DBt 
V ' \SdMx), xeD, v ' 



where the pair (<p,t/j) € H^, 1 ^ 2 (dD) x H 1 / 2 (9£>) is the solutions to 
'S D [iP](x)-S D [<p}(x)=h(x) 



8S d [tP] . . dS° D [tp] 



dv 



dh 



( x ) W = zrW 



for x e dD. (2.17) 



dv 
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Even if k = [i = or /i = oo, we have a similar representation: 

u K (x) = h(x)+S%[<p K }(x), xeD c ,k = 0,w. (2.18) 
When k = (i = 0, ipo satisfies 

-I + (KP D y)[ l p ] = -— on 3D, (2.19) 



and if fi = oo, then ipoo satisfies 



-I + (n)*) fox,] = on &D. (2.20) 

We emphasize that ip K e i/^ (dD). See, for example, [5] for details of the above mentioned 
representation of the solutions. 

A similar representation formula holds for the solutions to the boundary value problems. 
Let ft be a bounded Lipschitz domain in M. d containing D, which is also Lipschitz. Let u be 
the solution to 

V-CV s u = inQ, (2.21) 

with either the Dirichlet boundary condition u = f or the Neumann boundary condition 
^=gon dCl. Let 

h(x):=-[ T°(x-y)p^ (y)da(y)+f dr ° {X ~ v) u(y)da(y), x G ft. (2.22) 
Jan wo - Jan ov (y) 



I en 

Then u is represented as 



h(x)+S° D [<p}(x), xeft\D, 



u(x) = i v ' ' n x ' (2.23) 

\SnMx), xeD, 

where the pair (tp,ip) e H~ 1,2 (dD) x H~ 1 / 2 (dD) is the solutions to (j2T7t . 



3 Energy estimates 



Let 



J[u] := i f CVuiV'u+J / C°V s (u-h) : V s (u-h). 



(3-1) 



Here and throughout this paper A : B = . =1 aijbij for A = (a,j) and B = (fry). For the 
solution u to (|2.6|) . we prove that J[u] is bounded regardless of k and //. More precisely we 
prove the following lemma. 



Lemma 3.1 Let u be the solution to t2.6\) . If X < A for some constant A, then there is a 
constant C depending on A, but otherwise independent of ji and k, such that 



J[u] < C 



dh 

dv 



H- 1 / 2 (dD) 



(3.2) 



G 



As a consequence of Lemma 13.11 we have 



Lemma 3.2 Let tp be the potential defined in H2.16}) . If A < A for some constant A, then 
there is a constant C depending on A, but otherwise independent of fi and k, such that 



IMIff-i/2(3D) < C 



Oh 



H- 1 / 2 (dD) 



(3.3) 



Proof. Let v := u h. Then (|2 . 16|) yields v(x) = Sjj[(p](x) for x G Z? c . Thus, we have 
from ([2TT2j) 

Since \l + {K,° D )* is invertible on H^ 1 ^ 2 {dD), we have 



IMIff-V3(9£>) < C 



<9v 






+ 



H-V 2 (dD) 



Let ?/ be a function in H x / 2 {dD) satisfying J dD ri = and let w be the solution to 
Aw = in D c with w(x) = 0(|a;| 1_£i ) and w = r\ on 91?, so that the following estimate 
holds: 

||V s w|| L 2 (£) c ) < C||r?|| H i/2 (aD) . 

Since 



i0r7S 



V s v : V s wdx = 



dv 

,,d dvo 



r)da(x), 



(3.4) 



we have 



r dv 


■ r]da(x) 




JdD 9v 


< 







< C\\V s v\\ L 2 (D c ) \\r]\\ H i / 2 {dD) . 
Since n € iJ 1 / 2 (<9Z3) is arbitrary, we have from (|3.2[) 

ah 



dv 






+ 



H- 1 / 2 (dD) 



□ 



and so follows (I3.3[) . 

Proof of Lemma \3.1\ Let v := u h. It is known (see for example [16]) that v is the 
mmimizcr m W(R d ) of the functional 



r M -= \ I C(V s v + l D GV s h) : (V s v + l D GV s h), 



(3.5) 
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where G = I4 — (C 1 ) 1 C° and I4 is the identity 4-tensor. Note that 

= I I (C x V*(v + h) - C°V s h) : (V s (v + h) - (C x )- x C°V s h) 
2 Jd 



id 
1 

2 



;°V s v : V s v 



- [ C x V s v:V s v + f (C 1 - C°)V s v : V s h 
2 Jd Jd 

+ \ ( (C 1 -2C° + C°(C 1 )- 1 C°)V s h: V s h+i / C°V s v : V s v. (3.6) 
2 Jn 2 J D a 



Let 



(3.7) 

Then Voo £ W(R d ), and (voo + h)\ D e * which implies that V s (voo + h) = in D. So, we 
have from the first line in (13.61) that 



/[vco] = I f C°(C x )- x C°V s h : V s h + ~ / C°V s Voo : V s Voo . (3.8) 
2 Jd ^ Jd<= 



We then have 



J[u] = J[v + h] 



1 / C 1 V s v:V s v+ / C x V s v:V s h+^ f C x V s h : V s h 

2 Jd Jd 2 Jd 



:°V s v : V s v 



ID" 

= /[v]+ f C°V s v:V s h+ f C°V s h : V s h - i f C°(C x )- x C°V s h : V s h 
Jd Jd 2 J D 

= J[v]+ / C°V s u:V s h-i f C°(C x )- x C°V s h : V s h. 
Jd 2 7 D 

Since I[v] < /[voo], it follows from (|3.8p that 

J[u]</[ Voo ]+ f C°V s u : V s h - I [ C°(C 1 )- 1 C°V s h : V s h 
Jd 2 J D 



2 



D' 



D 

°V s Voo : V s Voo + / C°V s u : V s h. (3.9) 



Note that since A is bounded, we have 



f C°V s u : V s h = C( f \V s u\ 2 + f |V s h| 2 ) 
Jd Jd Jd 

<C( - f C x V 5 u:V s u+ f |V s h| 2 
Jd Jd 

< C fij[u] + / |V s h| ; 
Jd 



So, if ji is sufficiently large, then we have from (13.91) that 



J[u]<c(||V s v 0O ||i 2(Dc) + ||V s h|| 2 L2(I); 



for some constant C . Since 

||V S V co || i 2( £1 c) < C||<A>o||jf-i/2(d.D) < C 



dh 



dv 



H- 1 / 2 (dD) 



we have (|3.2j) when /i is large. 

When k and /i are bounded, we need a function which plays the role of Voo in the above. 
For that we use (po in (|2.18j) : define 



(3.10) 



It is worth emphasizing that vo is defined not only on D° but on M. d . Then one can show 
as above that 

J[u]<J[v ] + / C°V s v:V s h+ / C°V s h:V s h-i f C°(C 1 )- 1 C°V s h : V s h. 

J D J D 2 J D 



Using (|3.6p . one can see that 



Vv : V s v 



D 



I 1 - C°)V s v : V s h + - / C^Vh : V s h 



D 



\ [ C°V s v : V s v + f C°V s v : V s h. 

2 JD a J D 



(3.11) 



Since ^i±^ 



on dB, we have 



C°V s v : V s h = 



dh 



dv 
,id 9va 



:°V s v : V s v 



< 



c 



V s v : V s v + Ce I C°V s v : V s v 



for a (small) constant e. If e is sufficiently small, then we obtain by combining this with 

J[u] < C(||V s v || 2 L2(Rd) + ||V s h||i 2(D) ) (3.12) 
for some constant C independent of n and /i. Since 

IIWolU^co) < C\\^\\ H -i/2 {aD) < C 
we have (|3.2p . This completes the proof. □ 



dh 



dv 



H- 1 / 2 (dD) 
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4 Convergence of potentials and solutions 



Lemma 13.21 shows that the potential defined in (|2.16l) is uniformly bounded with respect to 
fi and A as long as A is bounded. We now prove the following lemma. 



Lemma 4.1 Let ip, ipo and ipoo be potentials defined by h2.16\) , \2.19\) and A2.20\) , respec- 
tively. 



(i) Suppose that A < A for some constant A. There are constants fix and C such that 



||V- ^oo||#-i/2(aD) 

for all fi > fii. 
(ii) There are constants 5 and C such that 



< 



C_ 



dh 



H- 1 / 2 (dD) 



(4.1) 



\\<P - ^0||h-V2( 9O ) < C(k + £i) 1/4 



dh 



dv 



(4.2) 



for all k, fi < 5. 



Proof. We may assume that || J^7-||.fr-i/2(aD) = 1. Let w = u so that w satisfies 



V ■ C°V s w = in D c , 
w - u|;q G W in D, 
r{x) = 0{\x\ l - d ) as|a;|^oo. 

Since V s (w — u) =0 in D, we have from Lemma T3. II that 

V s w\ 2 dx < — [ C x V s w : \7 s wdx 

D tyjD 



(4.3) 



1 



2/t 



— / ^V'u : X7 s udx < — . 



C 



(4.4) 



Let ipj, j = 1, . . . , d(d + l)/2, be a basis of "J as before, and let 



d(d+l)/2 



where f3j are chosen so that 



D 



(w-e)-Vi+V(w-e): =0, j = 1, . . . , d(d+ l)/2 



We then apply Korn's inequality (see for example JT9 1 ) to w — e to have 



|w-e|^ + |V(w-e)|^ <C |V s (w - e)\* < C / |V s w|- 



(4.5) 



(4.6) 
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for some constant C independent of It then follows from (|4.4[) that 

C 

l|w-e|| ff i (£)) < — , 

and from the trace theorem that 

C 

||w - e|| H i/a (aD ) < — 

for some constant C independent of /j,. By the strong convexity 
constant C such that 



(4.7) 
(4.8) 

of C°, there is a 



|V s w|| 2 L2(Z5c 



< C [ C°V s w : \7 s wdx 








■ w da(x) 


= -C [ — 

JdD d v Q 




• (w — e) da(x 

+ 


= -C [ — 
JdD dva 





dw 






e \\H 1 / 2 (dD) 




+ 


H- 1 / 2 {dD)'' 



where the second equality holds because of the orthogonality property (|2.10p . It then follows 
from P~5|l that 

C 



s„,||2 



|V s w 



We then obtain using the H^ 1 / 2 -!! 1 / 2 duality, divergence theorem on D c , Korn's inequality, 
and the trace theorem that 



and from (|4.9p that 



dw 




2 






+ 


H- 


-1/2 




dw 






9vq 


+ 



< C\\V s w\\ 2 L2(Dc)l 



< 



c 



Using the representations (|2.16|) and ()2.18|) . we have 

w(x) = u(x) - u oc (a;) = &b[tp - (f^x), x G D c 

Thus, ([2~T2|) yields 



dw 

dv a 



J + (£$,)* )[v -¥»«>] ondD. 



(4.10) 

(4.11) 
(4.12) 



So, gH]) follows from (fCTO]) . 

To prove (|4.2p . let v be as defined in (|3.10[) and let z := v — v in M. d . Then z = u u 
in D c and the following holds 



\V s z\ 2 dx < C 



•<0v7S 



V s z : V s zda; 



-C 



dz, 



= -C 



/" <9u 


/" du 




/an 


+ JdD £>vo 


+ -u ) 
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for some constant C > 0. Since 
du 



OD 



du Q 



du 



■ u = / C'V s u : V s u > 0, 

JD 



we have 



W s z\ 2 dx < C 



C 



du 

,,d d v o 
du 



c)D 



du 
°\ P 

IdD OV 



■ u 

• (v + h) 
(v + h) 



C / C'V'u: V s (v + h). 

JD 



By Cauchy's inequality, we obtain that 
J C x V s u: V s (v + h) < : V s u) 1/2 (y C^^vq + h) : V s (v + h)) 



1/2 



Thus, from (l3~2l it follows that 



y \V s z\ 2 dx <c(J C^^vo + h) : V s (v + h; 



1/2 



< CV^+7i||V s (v + h)\\ L 2 {D) < Cy/^+Jl 
for a constant C independent of k. Therefore, we arrive at 



dz 




dvo 


+ 



H- 1 / 2 (dD) 



< C||V s z|| L2(Dc) <C( K + fi)^ 4 . 



(4.13) 



Note that z = u — Uo = 5^ [p — ipo] in Z? c . So by the same reasoning as above we have (|4.2[) , 
and the proof is complete. □ 

As a consequence of Lemma |4. 11 we obtain the first main result of this paper. 

Theorem 4.2 Suppose that i2.2\) and 12. 3\) hold. Let u, Uoo and Uq be the solutions to 
[Kb)) . (EM and $7fy , respectively. 

(i) Suppose that A < A for some constant A. There are constants ji\ and C such that 



||u - Uoo||w-(H d ) < 

for all fjL> H\. 
(ii) There are constants 5 and C such that 



C 



dh 



dv 



H- 1 / 2 (dD) 



(4.14) 



\\u-U Q \\ w{Dn <C(K + ri 1/4 



dh 



dv 



H- 1 / 2 (dD) 



(4.15) 



for all k, fi < 5 . 
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It is not clear if the convergence rate, \x 1 / 2 and (k + /i) 1 ^ 4 , are optimal or not. 

Proof of Theorem \4.2[ Assume that \\-§^\\H- 1 / 2 {dD) = 1- Since u — u = Sjj[ ip — <po] on D c , 
dHHl) follows from ((Ml). 

Since u — = S^lip — tp^] on D c , we have 

C 

||u- Uoollw^c) < — . (4.16) 

Moreover, we have 

C 

||u - Uqo || #1/2(90) = W^dIV ~~ ^oo]||#i/2(aD) < C\\ip - <Poo||#-i/2(9D) < 

and hence 

llu-u^Uffi^) < — . (4.17) 
This completes the proof. □ 



5 Boundary value problems 

We now show that the results on the boundedness of the energy functional and on the 
convergence of solutions similar to the previous ones hold for the boundary value problems. 

Let fi be a bounded domain in R d and let D be an open subset in O. We assume that 
VL and D have Lipschitz boundaries and satisfy 



for some c > 0. We consider 



dist(£>, dQ) > c 



( V • CV s u = in n, 

— = g on dil, 

dv 

[u\ da e L%(dty, 



(5-1) 



where C is given by 



l^ + l^C . 



The Dirichlet problem can be treated in the exactly same way. 

The relevant energy functional for the boundary value problem is 



JnM ~ \ I CV s u:V s u. 
^ Jn 



(5.2) 



(5.3) 



Then the solution u to (|5.ip is the minimizer of Jji over H 1 ^) with the given boundary 
condition. Let be the solution when /i = oo (A is bounded). Then we have 



Jh[u] < Jn[uoo] = - / CV s Uoo : V s Uoo . 
z Jn 
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Since u^Id e "J, we have C 1 V s u DO = in D, and so, 

Jo[u] < \ [ C°V s Uoo : V s Uco < C ||g|| H -v»(«i) . ( 5 - 4 ) 



m\D 

where C is independent of /x, A. 

Using (|5.4p one can show as before that 



C 

||u- Uoolln^n) < — (5.5) 

for all fi > fix when A is bounded, and 

||u-uo|| ff i(n\D) <C(^ + M) 1/4 (5-6) 

for all k, fi < 5. Here Uo is the solution when k = \i = 0. We also note that (|5.4p together 

with Korn's inequality implies that ||u||^i/2(g n ) < C independently of fi, A. It then follows 
from flH22j) that 

\\Hm(n)<C (5.7) 

independently of n, A. 



6 Uniformity of asymptotic expansions 

Let us first recall the notion of elastic moment tensors (EMTs) associated to the inclusion D 
with the Lame constants ((m, A) when the background Lame constants are (Xq, /iq). Let a,ft£ 
N d be multi-indices, and let {efc}J? =1 be the standard basis of R d . For a £ N d , j = 1, . . . , d, 
let (ip{,ip J a ) be the solution to ([2T7) with h replaced by x a e r Here x a = x* 1 ■ ■ ■ x^ d . The 
EMT associated with D is defined by 

K P {D) = f x^ a da (6.1) 

JdD 

for a, (3 E N d and j — 1, . . . ,d. It is worth mentioning that is a vector, i.e., M° a g = 

( TO a/3 > TO a/3 > TO f/3 ) where 

m*(D) = / x^.^dtr. (6.2) 

If }ck) = = 1, we may write as m,j Pg for i, j,p, q = 1, . . . , d. It is known that (mij pq ) 
is an (anisotropic) elasticity tensor. See [12]. We emphasize that 

Wo\\H-V*{dD) < C (6-3) 

for some C independent of [i and k as long as A is bounded. 
Suppose that D is diametrically small and it is given by 

D = z + eB, (6.4) 



14 



where e is a small parameter representing the diameter of D, B is a reference domain 
containing 0, and Zq represents the location of D. Let u be the solution to (|5.ip and U be 
the background solution, i.e., the solution to 



V ■ C°V S U = 
dU _ 

LU| 8n eL|(3Q). 



in f2, 
on dVt, 



(6.5) 



Let N be the Neumann function which is the solution to 

'V ■ C°V s N{x,y) = -5 y (x)I 
0N(-,y) _ 1 



in O, 
on dfl. 



(6.6) 



,N( -,y) G for each y G ft, 



where I is the identity matrix. It is proved in [T^] (see also [7]) that the following asymptotic 
expansion holds on dft: for x G dft 

d d rf+l-|a| | Q | + | /3 | + d _ 2 

u(x)=U(x)-X)E E i«i ^[/ i ( 2o )^N( :E ^o)M^( J B)+0(e 2d ), 

3=1 |a|=l |/3|=1 a ' P ' 

where U = (C7"i, . . . , Uj). Our goal in this section is to show that this asymptotic formula 
holds uniformly in A and /x. More precisely we have the following result. 



u(«0 = U(x) - £ £ £ 

i=i M=i 101=1 



d d d+l-|a| e | Q | + |,g| +( i_ 2 



Theorem 6.1 Suppose that i2.2\) and \2.3\) hold. We have for x G dfl 

i 

, -9 a C/ J (^ )afN( : r,^)M^( J B)+E( : r), (6.7) 
|=i a ' p - 

where the error term satisfies 

||E|U~ (8n) < Ce 2d (6.8) 
/or some constant C independent of [i and A as Zong as A < A /or some constant A. 

We emphasize that (|6.7[) contains not only the leading order (e d ) term but also higher 
order terms up to e 2d_1 . The terms higher than e 2d are expressed in terms of not only EMTs 
but also interactions between the boundary <9f2 and the inclusion(s), and become much more 
complicated. 

Proof of Theorem 16.11 We closely follow the proof in [12] . By (|2 . 23[) , the solution u can be 
written in the form: 



u(x) 



S D [ip](x), 



xed\D, 

x G D, 



(6.9) 



where h is the function given by (|2.22j) and tp,ip G H 1 / 2 (dD) are the solutions to 
'S D [iP}(x)-S D [<p}(x)=h(x) 



dS D [i>] 



dp 



(x) 



.«-£<•> 



for x G dD. 



(6.10) 
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Let (p(x) := e(ys(zo + ex) and 4>{x) := eip(zo + ex). By a change of variables, (|6.10[) can 
be scaled into, for x £ dB, 



SbW\{x) - S%m(x) = h(z + ex) - 6, 



(12- 



Tn< 
2^T 



riLS 



dS B $] 



dv 



(x) 



dS B [0 



(6.11) 



dv 



(x) = tt- h ( z o +ex), 

OVq 



where 5d2 is the Kronecker delta function. By the Taylor expansion of h = (hi, . . . , hd), we 
have 



d d 



h(z + ex) = J2 J2 —d a h J (z )x a e J + 0(e d+1 ). 



(6.12) 



3=1 |a|=0 



Here the error term is independent of A, /i, and d a h(zo) is bounded independently of A, /i 
as long as A is bounded because of (|5.7|) . By (|6.12|) and the linearity of (|6.11|) we have 



d d 



^) = EE — 3 Q MzoV Q + o( e d+1 ) 



(6.13) 



3=1 |a| = l 



where ^ is given in the definition of EMT. Here we see from Lemma 13.21 that the error 
term in (|6.13l) is uniform with respect to A, as long as A is bounded. 
It is known that 



u(x) = U(x) - / N(x,y)ip(y)da(y), x € dQ 

JdD 



(6.14) 



(see [H]). Using (|6 . 1 3|) and the Taylor expansion of N(cc, y), we have 

oo 

N(as, z + ey)=Y, ^ W d^N(x, z Q )y^ x e dQ, 

and so, for x € 317, 

«(x)=U(a?)- / N(aj, z + ey)<p(y)e d - 2 da{y) 

JdB 



d d d+l-\u\ 
3=1 |a| = l 1,31 = 1 



*| + |0|+d-2 



[;~J! 



9«^(zo)^N(x,z )M^ + 0(e 2rf ). (6.15) 



The formula f|6. 15[) implies in particular that 

l|u-U|| ffl/w =0(e d ), 

where 0(e d ) is uniform with respect to A and /i. Since 



h(x) 



dT l {X , , V) "(y)Mv)~ I T°(x-y) S (y)da(y) 



XJ(x) 



dT°(x-y) 
dva{y) 



(u-V)(y)da(v), 
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we have 

d a h(z ) = d a XJ(z ) + O(e d ), (6.16) 
independently of A and /i. By substituting this into (|6.15p . we have (|6.ip . □ 



7 The case of multiple inclusions 

So far we deal with the case where the inclusion D is a simply connected inclusion, but 
the methods and results of this paper work even when there are multiple simply connected 
inclusions. 

Let us make a brief remark on the case when D has n disjoint simply connected compo- 
nents, say Di,... ,D n . In this case, The solution u to (|2.6[) is represented as 

{n 
3 = 1 (7.1) 
Sj[il>M](x), xeDj, j = l,...,n, 

where and Sj denote single layer potentials on dDj, and tpV) , ipW are the solutions to 



3 = 1 



8Sj[^ 



W-E 



9h ( 



for x e dDj, j = l,...,n. (7.2) 



If k is either or oo, we have a similar representation: 



u K (x) = h(x) + J2 S j[ l P% ) }( x )> xeD c , K = 0,OO, 

3 = 1 



(7.3) 



where ip)z satisfies appropriate integral equations. We can show in a similar way that 



eh^-^iih-v-^j^-Ie 

3 = 1 V» J = 1 



dh 



and 



dh 



(7.4) 



8vq 



(7.5) 



E H^ (J) " WWs-V'vn,) < C(k + ^) 1/4 E 

So, Theorem 14.21 and Theorem 16. II are valid even if D has several components. 

It is worth mentioning that if there are multiple inclusions, the convergence depends on 
the distance between inclusions since |Vu| may be arbitrarily large as the distance between 
inclusions tends to zero. This fact was proved in, for example, |11U301IT5] for the conductivity 
problem. (See [3] for an extensive list of recent papers on this problem.) For the elasticity 
problem, it is shown in |24j by numerical computations that |Vu| may blow up as the 
distance between inclusions tends to zero. 
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8 Conclusion 

In this paper we have used layer potential techniques to prove uniform convergence in an 
appropriate function spaces of solutions to the Lame system as the bulk and shear moduli 
tend to extreme values (zero or oo) provided that the compressional modulus is bounded. 
Making use of this result, we have shown that the asymptotic expansion of the solution due 
to the presence of diametrically small inclusions is uniform with respect to the bulk and shear 
moduli. These results are obtained under the assumption that the Lame parameters of the 
background and inclusions are constant. We expect that the same results hold even if the 
background Lame parameters are not constants, but variables, even though the methods of 
this paper do not apply to that case. Another interesting case is when the inclusion is thin, 
and the thickness tends to zero |15j . In this case we expect that the asymptotic expansion 
is not uniform. 
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